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^ ' We show that QCD Dirac spectra well below Aqcdi both at zero and at nonzero 

QP ' chemical potential, can be obtained from a chiral Lagrangian. At nonzero chem- 

ical potential Goldstone bosons with nonzero baryon number condense beyond a 
critical value. Such superfluid phase transition is likely to occur in any system 
with a chemical potential with the quantum numbers of the Goldstone bosons. We 
discuss the phase diagram for one such system, QCD with two colors, and show 
the existence of a tricritical point in an effective potential approach. 
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i-S^ ■ 1. Introduction 

. ^ For strongly interacting quantum field theories such as QCD a complete 

■ nonperturbative analysis from first principles is only possible by means of 

^ I large scale Monte Carlo simulations. Therefore, partial analytical results in 

some parameter domain of the theory are extremely valuable, not only to 
provide additional insight in the numerical calculations, but also as an in- 
dependent check of their reliability. This has been our main motivation for 
analyzing such domains. The principle idea we have been pursuing is based 
on chiral perturbation theory [ 1,2]: because of confinement and the spon- 
taneous breaking of chiral symmetry, the low-energy chiral limit of QCD 
is a theory of weakly interacting Goldstone bosons which are described by 
a chiral Lagrangian that is completely determined by the symmetries of 
QCD. This idea can be applied to the QCD Dirac spectrum which can be 
extracted from the valence quark mass dependence of the chiral condensate. 
The valence quark mass is not a physical parameter of the QCD partition 
function and can be chosen in a domain where the valence quark mass de- 
pendence of the QCD partition function can be described to an arbitrary 
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accuracy by a corresponding chiral Lagrangian. If the Compton wavelength 
of Goldstone bosons containing only valence quark masses is much larger 
than the size of the box the low-energy effective theory simphfies even much 
further [ 3] . Then only the zero momentum component of the Goldstone 
fields has to be taken into account so that the valence quark mass depen- 
dence of the QCD partition function is given by a unitary matrix integral. 
This idea was first applied to the QCD partition function [ 4] with quark 
masses of order m ~ 1/V{tp^p) (with {■iptp) the chiral condensate and V the 
volume of space-time). However, we emphasize that for physical values of 
the quark masses and volumes, a part of the Dirac spectrum, as probed 
by the valence quark mass, is always in this mcsoscopic domain of QCD. 
More precisely, using the Gell-Mann-Oakes-Renner relation (with F the 
pion decay constant), in the domain 

and T/V4A^cD»l, (D 

the kinetic term in the chiral Lagrangian can be ignored and the valence 
quark mass dependence of the QCD partition function reduces to a unitary 
matrix integral [ 3,5]. This integral is equivalent to a chiral Random Matrix 
Theory in the limit of large matrices [6,7]. The second condition ensures 
that excitations of the order Aqcd decouple from the low-energy sector of 
the partition function. 

At nonzero baryon chemical potential the Dirac spectrum is scattered 
in the complex plane. However, at a sufficiently small nonzero baryon 
chemical potential and finite physical quark masses, the Dirac spectrum in 
the phase of broken chiral symmetry is still described by a partition function 
of Goldstone bosons containing valence quarks [ 8] . In order to eliminate 
the fermion determinant containing the valence quarks, one has to calculate 
the valence quark mass dependence of the chiral condensate in the limit of 
a vanishing number of valence quarks. The existence of this limit requires 
the introduction of conjugate antiquarks [ 9,10], resulting in the appearance 
of Goldstone bosons with nonzero baryon number containing only valence 
quarks. They condense if the chemical potential exceeds their mass. In 
terms of the Dirac spectrum this phase transition is visible as a sharp 
boundary of the locus of the eigenvalues. 

Such phase transition to a Bose condensed phase is likely to occur in any 
theory with a chemical potential with the quantum numbers of Goldstone 
bosons. For example, for QCD with two fundamental colors [ 11,12] or for 
adjoint QCD with two or more colors [ 12], the lightest baryon is a Gold- 
stone boson. A transition to a Bose condensed phase occurs for a chemical 
potential larger than the mass of this boson. Other examples are pion con- 



3 



densation, which may occur for a nonzero isospin chemical potential [ 13], 
and kaon condensation which may occur for a nonzero strangeness chemical 
potential [ 14] . If the mass of the Goldstone bosons and the chemical po- 
tential are both well below Aqcd, such phase transition can be described in 
terms of a chiral Lagrangian. We have analyzed such Lagrangian for QCD 
with two colors at nonzero temperature and chemical potential [ 12,15,16]. 
In an effective potential approach we have found a tricritical point [ 16] in 
agreement with recent lattice QCD simulations [ 17]. 

We start this lecture by discussing QCD Dirac spectra at zero chemical 
potential and explaining its description in terms of a chiral Lagrangian. 
In section 3 we analyze QCD Dirac spectra at nonzero chemical potential. 
The phase diagram of QCD with two colors at nonzero temperature and 
chemical potential is discussed in section 4 and concluding remarks are 
made in section 5. 



2. Dirac Spectrum at Zero Chemical Potential 

The Euclidean QCD Dirac operator is givc;n by 

iD = -f^{d^ + iA^), (2) 

where the 7^ are the Euclidean gamma matrices and the arc SU{Nc) 
valued gauge fields. The Dirac spectrum for a fixed gauge field configuration 
is obtained by solving the eigenvalue equation 

iD(j)k = iXk<i>k- (3) 

In a regularization scheme with a finite number of N eigenvalues, the av- 
erage spectral density is defined by 

N 

p(A) = (^(5(A-Afe))QCD, (4) 
fe=i 

where the average (• • •)qcd is over gauge field configurations weighted by 
the Euclidean QCD action. As a result of the averaging we expect that p(A) 
will be a smooth function of A. Because of the involutive automorphism 
"fsiD'y^ = —iD the Dirac operator can always be represented in block-form 
as 

If ly is a square matrix the nonzero eigenvalues of iD occur in pairs ±Afe. 
For nonzero topological charge the total number of zero eigenvalues is given 
by the diflference of the the number of right-handed modes and left-handed 
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modes. In that case, the matrix is a rectangular matrix with the absolute 
value of the difference between the number of rows and columns equal to the 

topological charge. For very large values of A the Dirac spectrum converges 
to the free Dirac spectrum so that the spectral density given by p(A) ^ VX^. 
The smallest nonzero eigenvalue, Amin, is of the order of the average level 
spacing and is thus given by 

A..„ = AA = -^. (6) 



2.1. Spontaneous Chiral Symmetry Breaking and 
Eigenvalue Correlations 

The chiral condensate is given by 



(tp^p) = lim lim lim tt ( Tr- 



A— >oo rn— >0 V— >oo V \ iD 



m 



QCD 



1 2mp{\) _ 
= hm lim lim — / ^ ' . (7) 

A-^oo m^O y^oc V 

The limit m ^ is taken before A ^ oo to eliminate divergent contribu- 
tions from the ultraviolet part of the Dirac spectrum (the ultraviolet cutoff, 
A, may also appear in the spectral density). Because of spontaneous break- 
ing of chiral symmetry, the limit V ^ oo cannot be interchanged with the 
limit TO ^ in (7). If the chiral condensate is nonzero the limits to — > 0+ 
and to ^ 0^ have opposite signs. This can only happen if p(0) V. If we 
expand the spectral density as 

p(A)=p(0+)+ai|A|+a2A2 + ..., (8) 

we obtain Banks-Casher formula [ 18] 

(ViV)= lim (9) 

In this article we avoid taking limits by mainly focusing on finite values of 
TO, V and A. 

Let us now consider the QCD partition function Z{mf), 

Zimf) = ([[Y[{iXk+mf))YM, (10) 

/ k 

where (• • •)ym denotes averaging with respect to the Yang-Mills action. 
Because in the thermodynamic limit the derivative of the partition function 
with respect to to/ has a discontinuity across the imaginary axis, we expect 
that its zeros are also located on the imaginary axis as well and, for finite 
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volume, are spaced as 1/1^. This average can also be written as an average 
over the joint eigenvalue distribution 

p(Ai, A2, • • •) = (-^(Ai - Xt)5{\2 -Xi)-- •)qcd, (11) 

where are the eigenvalues of the Dirac operator for a given gauge field 
configuration A. This results in 

Z{mf)= f p{Xi,X2,---)'[l]l{iXk+m). (12) 

f k 

If the eigenvalues are uncorrelated the joint eigenvalue distribution factor- 
izes into one-particle distributions 

p(Ai,A2,---) = Pi(Ai)pi(A2)--- (13) 

and the partition function is the product of TV identical factors. For exam- 
ple, for Nf = 1, in the sector of zero topological charge, we obtain 

Z(m) = ((A2)i+mT, (14) 

where (•••)! is the average with respect to the one particle distribution 
(which in this case is the average spectral density of the QCD Dirac op- 
erator). Therefore Z(m) is a smooth function as m crosses the imaginary 
axis along the real axis and chiral symmetry is not broken. We conclude 
that the absence of eigenvalue correlations implies that chiral symmetry is 
not spontaneously broken, or conversely, if chiral symmetry is broken spon- 
taneously the eigenvalues of the Dirac operator are necessarily correlated. 
The question we wish to answer is what are these correlations. 

2.2. Low Energy Limit of QCD 

Because of confinement the chiral limit of QCD at low energy is a theory of 
weakly interacting Goldstonc bosons. For small values of the quark masses 
m/ and chemical potentials /z/ the QCD partition function coincides with 
a partition function of Goldstone bosons: 

ZQCD{mf,iJ.f,9) ~ ZGoid{mf,iJ,f,9), (15) 

where 6 is the vacuum ^-angle. Up to phenomenological coupling constants, 
the mass dependence of ^GoM is completely determined by the symmetries 
and transformation properties of the QCD partition function. In particular, 
both partition functions have the same low mass expansion. Equating the 
coefficients of powers of the quark masses leads to sum-rules for the inverse 
Dirac eigenvalues [ 4] . To derive them we consider the Fourier components 
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of the 6 dependence which are just the partition function in a given sector 
of topological charge, 

Zu{---) = ^ fj dee'^'Z{---,e). (16) 

As an example, let us consider the case Nf = 1, iif =0 and u = Q. In this 
case there arc no Goldstone bosons and the the mass dependence of the 
partition function for ^ = is given by 

The 6 dependence id obtained from the substitution m — > me'^. For the 
sector of zero topological charge we thus find 

27r 

mVE cos 



1 r^^^ 

{{X^+m^)),=o^ — dee 



l + imV^E2 + .... (18) 



This result in the sum rule [ 4] 



In fact, an infinite number of sum rules can be derived for the partition 
function of QCD and QCD-like theories with spontaneous chiral symmetry 
breaking [ 4,19,20,21,22]. Nevertheless, these sum rules are not sufficient to 
determine the Dirac spectrum. 



2.3. Resolvent 

In order to derive the QCD Dirac spectrum we introduce the resolvent 

Here, z is a complex 'valence quark mass' which does not occur inside the 
fermion determinant that is included in the average. The spectral density is 
obtained from the discontinuity of the resolvent across the imaginary axis, 

p{X) = ^{G{^\ + e)-G{^X-e)) 

= ^{G{iX + e)+G{-i\ + e)). (21) 
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The resolvent can be obtained [ 23,24,25,26,27,28,29,30] from the generating 
function Z^pgctiz, z',mf), 

G{Z) = Z,peet(^, z', mf)l,^^ , (22) 



with 

Zspect{z,z',mf) 



\ / ' YM 



The variable 2; is a parameter that probes the Dirac spectrum and can be 
chosen arbitrary small. For z, z', nif ^ Aqcd this partition function can 
be approximated arbitrarily well by a chiral Lagrangian which is completely 
determined by the symmetries of the QCD partition function. In addition to 
fermionic quarks, this partition function also contains bosonic ghost quarks. 
The corresponding chiral Lagrangian therefore includes both bosonic and 
fermionic Goldstone bosons with masses given by 2Re{z)T, / , Re{z + 
z')T,/F^, Re{z + mf)'S/F'^, etc., as given by the usual Gell-Mann-Oakes- 
Renner relation. 

The inverse fermion determinant can be written as a convergent bosonic 
integral provided that Ke{z) > 0: 



det{iD + z) 



d(j>d(j>*e-'^'^'^+^^'^. (24) 



The convergence requirements restrict the possible symmetry transforma- 
tions of the partition function. For example the axial U{1) transformation 
which in the fermionic case is given by 

i>R ^ e-^ViJ, i>L ^ e+^Vi, (25) 

would violate the complex conjugation structure of the bosonic integral 
with (f)ii = and ^i, = (^\. Instead, the allowed Va^) transformation is 

4>R ^ e-^iJ, 4>L ^ e'cj>L, (26) 

with s a real parameter. The Goldstone manifold is therefore not given by 
the super- unitary group SU {Nf + 1 1 1) but rathca- by its complexified version 
that reflects the convergence requirements of the bosonic axial transforma- 
tions [ 29,30]. We will denote this manifold by SU{Nf + l\l) and an explicit 
parameterization for the simplest case, Nf = 0, will be given below. Vector 
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flavor symmetry transformations are consistent with the complex conjuga- 
tion properties of the bosonic integral. This symmetry group is thus given 
hy SU{Nf + l\l). 

In the chiral limit the mass dependence of generating function (23) can 
be obtained from a chiral Lagrangian determined by its symmetries and 
transformation properties. It is given by 

L = Stvd^Ud^U-^ - i(V;V')Str(M(C/ + U-^)), (27) 

and the corresponding partition function roads 

Z^= [ dU{x)Sdet''Uoe-^'^^''^. (28) 

Ju{Nf + l\l) 

Because v is the global topological charge only the zero momentum compo- 
nent of U, denoted by Uq, appears in the argument of the supcrdctcrminant. 
In the chiral limit, QCD is flavor symmetric so that the kinetic term of the 
chiral Lagrangian should be flavor symmetric as well. Therefore, the pion 
decay constant of the extended flavor symmetry is the same as in QCD. 
The mass matrix is given by M = diag(?7?,, • • • . to, z, z'). 

If z <C Too = F^/{%l>^)VV the fluctuations of the zero momentum modes 
are much larger that the fluctuations of the nonzero momentum modes, 
which then can be ignored in the calculation of the resolvent. More phys- 
ically, this condition means that the Compton wavelength of Goldstone 
bosons containing ghostquarks with mass z or z' is much larger than the 
size of the box. In condensed matter physics, the energy scale TOc is known 
as the Thouless energy and has been related to the inverse diffusion time 
of an electron through a disordered sample [31]. 

In the Dirac spectrum we therefore can distinguish three different energy 
scales, the smallest eigenvalue Amin, the Thouless energy to,, and the QCD 
scale Aqcd- On mass scales well below Aqcd the mass dependence of the 
QCD partition function is given by the chiral Lagrangian. For mass scales 
well below the Thouless energy only the zero momentum modes have to 
be taken into accoimt. However, for masses not much larger than Xmim a 
perturbative calculation breaks down and the group integrals have to be 
performed exactly. An interesting possibility is if Amin and rric coincide 
which may lead to critical statistics [ 32] . 

In the zero momentum limit, it is straightforward to calculate the inte- 
grals over the superunitary group. The simplest case is the quenched case 
{Nf = 0) where U can be parameterized as 
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with a and (3 are Grassmann variables, 9 € [0, 27r] and s G (—00,00). In 
terms of the rescaled variable u = zVl-tptp), one obtains the resolvent 

G{u) 



u{Ka{u)Ia{u) + Ka-l{u)Ia+l{u)) + -, 



(30) 



where a — Nf + \v\. Prom the definitions of the modified Bessel functions it 
is clear that the compact /noncompact parameterization of the supcrunitary 
group is essential. The microscopic spectral density is obtained from the 
discontinuity of the resolvent and is given by 



p((/v'(tn')) c 



where Q = XV{'tjjil)). 



^{J^iO - Ja+i(C)Ja-l(C)) + l^SiO, (31) 



1.2 
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Figure 1. The valence quark mass dependence of the chiral condensate T,{mv) plotted 
as T,(z/m,v)/'B versus rrivVT,. The dots and squares represent lattice results by the 
Columbia group [ 33] for values of /3 as indicated in the label of the figure. (Figure 
adapted from ref. [3]). 



2.4. Lattice Results 

The properties of the Dirac spectrum have been analyzed in many lat- 
tice QCD simulations [33,3,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48] 
[ 49,50,51,52,53,54] and have been found to be in complete agreement with 
the conclusions of the previous section. We only show three representative 
examples. 
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In Fig. 1 we show the valence mass dependence of the chiral condensate 
as calculated by the Columbia group [ 33] . In this figure the valence quark 
mass is denoted by niv and S = {'ip^p). Our variable u in (30) is thus 
given by rriyVE and S(mt,) should be identified with G{u). The reason 
that the lattice data agree with the quenched approximation is that the 
sea-quark masses in the lattice calculation are much larger than the valence 
masses. The topological charge is zero because the instanton zero modes are 
completely mixed with the nonzero modes due to the lattice discretization. 

Because the valence quark mass dependence agrees with (30) the corre- 
sponding lattice QCD microscopic spectral density should agree with (31). 
This was shown by two independent calculations [ 42,43]. In fig. 2 we show 
results for an S'' lattice with quenched staggered fermions [ 42] . 




5 10 15 20 

Figure 2. The microscopic spectral density for quenched QCD with three colors. The 
solid curve represents the analytical result (31) for a = 0. (Figure adapted from ref. [ 42]) 



In Fig. 3 we show the disconnected chiral susceptibility defined by 

= TV ( {i\k+m,){i\i+m,))~N iXk+m,) "^^^^ 

This quantity can be obtained from the two-point spectral correlation func- 
tion but can also be directly computed in chPT [ 55,56,29,54]. The dashed 
curve represents the result obtained from taking into account only the zero 
momentum modes whereas the solid curve is obtained from a perturbative 
one-loop calculation. Also in this figure S = {'tp'ip). This figure clearly 
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demonstrates the existence of a domain where a perturbative calculation 
can be applied to the zero momentum sector of the theory. 
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Figure ?>. The disconnected susceptibility for quenched SU(3) with staggered fermions 
(solid points). The solid curve represents the prediction from chPT, and the dashed one 

is the exact result for the zero momentum approximation to the chiral susceptibility. 
(Note the dashed line is hidden by the data points for u < 10.) (Figure taken from ref. 
[ 54].) 




2.5. Chiral Random Matrix Theory 

Correlations of Dirac eigenvalues on the scale of the average level spacing 
are completely determined by the zero mode part of the partition function 
which only includes the mass term and the topological term of the chiral 
Lagrangian. This raises the question of what is the most symmetric the- 
ory that can be reduced to this partition function. The answer is chiral 
Random Matrix Theory in the limit of large matrices. This theory is in- 
variant under and additional Uji(n) x UL^n + i/) group (with n x + i/) 
the size of the nonzero blocks of the Dirac matrix). Because of this much 
larger symmetry group, all correlation function of the eigenvalues can be 
obtained analytically, often in a much simpler way than by means of the 
supersymmetric generating functions for the resolvent. 

Before defining chiral Random Matrix Theories, we have to introduce 
the Dyson index of the Dirac operator. It is defined as the number of 
independent degrees of freedom per matrix element and is is determined by 
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the anti- unitary symmetries of the Dirac operator. They are of the form 

[AK,iD]=0, (33) 

with A unitary and K the complex conjugation operator. As shown by 
Dyson [ 57] , there are only three different possibilities within an irreducible 
subspace of the unitary symmetries 

i) There are no anti — unitary symmetries, 
ii) {AKf = l, 

Hi) {AKf = -l. (34) 

In the first case the Dirac operator is complex and the Dyson index is 
= 2. In the second case it is always possible to find a basis in which the 
Dirac matrix is real and the Dyson index is = 1. In the third case it is 
possible to express the matrix elements of the Dirac operator into selfdual 
quaternions and the Dyson index is (3n = 4. The first case applies to QCD 
with three or more colors in the fundamental representation. The second 
case is realized for QCD with two colors in the fundamental representation, 
and the third case applies to QCD with two or more colors in the adjoint 
representation. 

Chiral Random Matrix Theory is a Random Matrix Theory with the 
global symmetries of the QCD partition function. It is defined by the 
partition function [6,7] 

Z,{mi,---,ms)= / dW^ JJdet"(i£' + m/)e-*'^^('^^*), (35) 

/=i 

where the Random Matrix Theory Dirac operator is defined by 

and is an n X (n + i^) matrix so that iD has exactly v zero eigenvalues. In 
general the probability potential is a finite order polynomial. However, one 

can show [58,59,60,61,62,63] that correlations on the scale of the average 
level spacing do not depend on the details of this polynomial and the same 
results can be obtained much simpler from the Gaussian case. Depending 
on the Dyson index of the Dirac operator we have three different possibil- 
ities, the matrix elements of W arc real, complex or self-dual quaternion 
for = 1, 2, 4, respectively. The corresponding Gaussian ensembles are 
known as the chiral Gaussian Orthogonal Ensemble (chGOE), the chiral 
Gaussian Unitary Ensemble (chGUE) and the chiral Gaussian Symplectic 
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Ensemble (chGSE), in this order. Together with the Wigner-Dyson Ensem- 
bles and four ensembles that can be applied to superconducting systems, 
these ensembles can be classified according to the Cartan classification of 
large symmetric spaces [ 64] . 

The results (30) and (31) quoted in the previous section were obtained 
first by means of standard Random Matrix Theory methods [ 55,3]. 

3. Dirac Spectra at Nonzero Chemical Potential 

Quenched lattice QCD Dirac spectra at /i 7^ were first obtained numer- 
ically in the pioneering paper by Barbour et al. [ 65] and have since then 
been studied in several other works [ 66,67,68,69]. Since the Dirac opera- 
tor has no hcrmiticity properties at /i 7^ its spectrum is scattered in the 
complex plane. However, it was found [ 65] that for not too large values of 
the chemical potential the spectrum is distributed homogeneously inside an 
oval shape with a width proportional to fj,^. In this section we will explain 
these results in terms of a chiral Lagrangian for phase quenched QCD at 
nonzero chemical potential. 

3.1. Spectra of Nonhermitian Operators 

The spectral density of a nonhermitian operator is defined by 



Often it is useful to interpret the real and imaginary parts of the resolvent 
as the electric field in the plane at point z from charges located at Afe. 

Since the fermion determinant is invariant for multiplication of the Dirac 
operator by an unimodular matrix, one could analyze the spectrum of var- 
ious Dirac operators. The Dirac operator that is of interest is the one 
with eigenvalues that are related to an observable. For example, the Dirac 
operator in a chiral representation has the structure 




= -9.*G(z), 



(37) 



where the resolvent G{z) is defined by 




(38) 




(39) 
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In terms of its eigenvalues, the chiral condensate is given by 

k 



(40) 



If we are interested in the baryon number, on the other hand, we consider 
the Dirac operator 



iW m 
m iW^ 



(41) 



which satisfies the relation det(i_D^ + /i) = det(ii3 + m). In terms of its 
eigenvalues fik the baryon density is given by 



riB 



(42) 



QCD 



Finally, let us consider QCD at nonzero isospin chemical potential. In this 
case the fermion determinant is given by 



det 



m iW + Hi 
iW^ +1-11 m 



det 



m iW — Hi 



(43) 



which can be rewritten as the determinant of the antihermitian matrix 
/ -m iW - 

iW^ - HI -m 

m iW + Hi 
\iW^+Hi m / 

In terms of its eigenvalues ink , the pion condensate is given by 



(44) 



(tt) 



(45) 



QCD 



where is the source term for the pion condensate. 



3.2. Low Energy Limit of Phase Quenched QCD 

The generating function for the quenched Dirac spectrum is given by the 
replica limit {Nf 0) of phase quenched QCD partition function [ 10] 
defined by 

Z = {[det{iD + Z + M7o)det*(i£> + z + At7o)]^OQCD 

= {[det{iD + Z + M7o) det(i£> + z* - M7o)]'^Oqcd. (46) 

Since this is a partition function of quarks and conjugate anti-quarks we 
can have Goldstone bosons with nonzero baryon number. For a chemical 
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potential equal to half the pion mass we thus expect a phase transition 
to a Bose condensed phase. For a quark mass much less than Aqcd this 
phase transition can be described completely in terms a chiral Lagrangian. 
In nonhcrmitian Random Matrix Theory, the technique to determine the 
spectral density by analyzing a corresponding Hermitian ensemble is known 
as Hermitization [ 70,71], 

The chiral Lagrangian is again determined by the symmetries and the 
transformation properties of the QCD partition function. These can be 
made more explicit if we rewrite the fermion determinant as 

, / Ml d + BR\ 
^''[-d^ + B, M, j' 

where Mi = M2 = dia.g{z, - ■ ■ , z, z* ,■ ■ ■ , z*) and Bl = Br = 
diag(/x, • • • , /i, — /U, • • • , — /i). For z = n = our theory is invariant un- 
der SUlC^Nj) X SUR{2Nf). For 2; ^ and ^ ^ this invariance can be 
restored if the the mass and chemical potential matrices are transformed 
as [ 2,8,13] 

Ml ^ VrMiV^\ Br ^ VrBrV^\ (48) 
M2 ^ VlMiV^\ Bl ^ VlBlV^\ (49) 

However, since Br(^ arc a vector fields we can achieve local covariance by 
transforming them according to 

BL^VL{d„ + BL)V^\ 

BR-*VR{do + BR)V^\ (50) 

In the effective Lagrangian local covariance is obtained by replacing the 
derivatives in the kinetic term by a covariant derivative given by [ 2] 

a^s ^ v^s = a^s - s^s + t,Br. (51) 

This results in the chiral Lagrangian 

i= — TrV^SV^Et - -Tr(MiSt + M2S). (52) 

In our mean field analysis to be discussed below we only need the static 
part of this Lagrangian which is given by [ 8] 

j^stat ^ —n^Ti BrY,Bl^'i - -Tr(MiS^ + MzS). (53) 
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3.3. Mean Field Analysis 

In this subsection we describe the mean field analysis [ 8] of the static La- 
grangian (53). In phase quenched QCD, baryonic Goldstone modes contain 
a quark with mass z and a conjugate antiquark with mass z*. According 
to the GOR relation their mass is given by 

M^=^-i±^. (54) 

If the chemical potential is less than M/2 only the vacuum state contributes 
to the QCD partition function. This results in 

Z = e^(^+-*)G. (55) 

We then find the following result for the resolvent and the spectral density 

G{z) = G, p{z) = 0, for iJL < M/2. (56) 

For /u > M/2 the baryonic Goldstone modes condense resulting a non- 
trivial vacuum field which can be obtained from a mean field analysis. The 
mass term and the chemical potential term in the static Lagrangian are 
respectively minimized by* 

Sc= and S,= f'j'^V (57) 



,0 1/ VI 0, 

A natural ansatz for the minimum of the static Lagrangian (53) is thus 
given by 

E = Sc cos a + Srf sin a. (58) 

An effective potential for a is obtained by substituting this ansatz into the 
static Lagrangian. It is given by 

L{a) = fi^F^Nf (cos^ a - sin^ a) - GNf{z + z*) cosa. (59) 

This potential is minimized for a given by 

< ^''Ipf^ ■■ sina = 0, L{d) = -GN f{z + z*), 

Prom the free energy at the minimum we easily derive the resolvent and 
the spectral density (see Fig. 4 in units with 2/x^F^/G = 1) 

m'<%#^: G{z)=G, p(A) = 0, 



(60) 



*The minimum is not unique which leads to massless Goldstone bosons in the con- 
densed phase. 
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We conclude that the Dirac eigenvalues are distributed homogeneously in- 
side a strip with width ~ /z^ in agreement with the numerical simulations 
[ 65]. For a discussion of correlations of eigenvalues of a nonlicrniitian op- 
erator we refer to the specialized literature [ 74,75,70,77,76,78,79]. 
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Figure 4. The distribution of eigenvalues of the Dirac operator in the complex x-plane. 
The resolvent given by eq. (61) is represented by the dotted curve. 



In [ 13] this analysis was applied to the problem of QCD at finite isospin 
density with a partition function that coincides with the phase quenched 
QCD partition function (46) [ 72]. In that reference [ 13] it was also shown 
that the ansatz (58) is a true minimum of the static Lagrangiaii. However, 
although we believe that it is an absolute minimum, this has not yet been 
shown. 

4. Real QCD and Nonzero Chemical Potential 

The analysis of the previous section can be repeated for any theory with 
a chemical potential with the quantum numbers of the Goldstone bosons. 
Both for QCD with two colors in the fundamental representation and for 
QCD with two or more colors in the adjoint representation, a baryon has 
quark number two and is a boson. For broken chiral symmetry some of these 
baryonic states are Goldstone bosons so that Bose-Einstein condensation is 
likely to occur if the baryon chemical potential surpasses the mass of the 
Goldstone bosons. For QCD with three or more colors in the fundamental 
representation we expect a similar low energy behavior if we introduce a 
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chemical potential for isospin [72,13] or strangeness [ 14] leading to pion 
condensation or kaon condensation, respectively. Below we only discuss 
QCD with two colors. 

4.1. QCD with No = 2 

For simplicity, let us consider QCD with both two colors and two flavors. 
In that case diquark mesons appear as flavor singlet. We thus have five 
Goldstone bosons, three pions, a diquark and an anti-diquark^. Because 
SU{2) is pseudo-real, the flavor symmetry group is enlarged to 5'C/(4). 
The quark-antiquark condensate breaks this symmetry spontaneously to 
Sp{4:) [ 20,73,80]. We can again write down a chiral Lagrangian based on 
this symmetry group. Also in this case we find a competition between two 
condensates, and in the Bose condensed phase, the chiral condensate rotates 
into a diquark condensate for increasing values of the chemical potential as 
in (58). The mean field analysis proceeds in exactly the same way as in 
previous section. For the chiral condensate we obtain 

2 

_ _ 777,-^ 

li>m^/2, = (VV')o^. (62) 

In Fig. 5 we show that our predictions agree with lattice simulations by 
Hands et al. [81]. The simulations were done for a 4^ x 8 lattice for 
SU{2) in the adjoint representation and staggered fermions which is in the 
same symmetry class as QCD with two colors in the fundamental repre- 
sentation. A similar type of agreement was found by several other groups 
[ 82,83,84,85,86] 

Results for QCD with two colors in the fundamental representation 
obtained in [ 84] are shown in Fig. 6. Again we find good agreement with 
the mean field results (62). Furthermore, if we plot the same data versus 
TO^/4/i^ the curve reminds us of the resolvent for QCD in phased quenched 
QCD after transforming the z dependence of the resolvent at fixed /i into a fj, 
dependence at fixed z (see Fig. 4). Since the condensate can be interpreted 
as the electric field at the quark mass due to charges at the position of the 
eigenvalues we have no eigenvalues for /i > and for a narrow strip 

along the / 4/U^ axis. In the remaining region the Dirac eigenvalues are 
distributed homogeneously. The absence of eigenvalues close to the /4/u^ 



tPor QCD in the adjoint representation, the diquarks appear as flavor triplet. For two 
flavors this results in three pions, three diquarks and three anti-diquarks in agreement 
with spontaneous symmetry breaJfing according to 5(7(4) — > 0(4) 
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Figure 5. The chiral condensate versus n/(m-^ /2) for QCD with two colors in the adjoint 
representation (Figure taken from ref. [ 81]). 

axis is a signature[ 87] of /3_d = 4. Indeed, such behavior has been identified 
both numerically [ 87] and analytically [ 77] . 



4.2. Beyond Mean Field 

One of the recurring questions in the study of phase transitions is the 
stability of the mean field analysis. In the following, we carry out a next- to- 
leading order study of the second order phase transition found at the mean- 
field level. Additional details can be found in [15,16]. We will concentrate 
on the free energy of the Bose condensed phase close to the mean-field 
critical chemical potential jic = M/2, with the leading order pion mass 
given by the GOR relation: = GmJF'^. 

The chiral Lagrangian (52) contains the two operators that have the 
lowest dimension in momentum space and that are invariant under local 
flavor transformations. There are, of course, many operators of higher 
dimension that fulflU these symmetry constraints. One has to introduce a 
systematic power counting to account for their relative importance [ 2] . Our 
power-counting scheme is the same as the one used in chiral perturbation 
theory, extended to include the chemical potential: p ~ /x ~ M ~ ^Jfn^, 
where p is a Goldstone momentum. The leading order chiral Lagrangian 
contains all the operators of order that fulfill the symmetry constraints. 
The next-to-leading order chiral Lagrangian contains all suitable operators 
of order p^. In general, for any Nj, there are ten such operators that 
contribute to the free energy. They are made out of traces of VySVi/S^ 
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Figure 6. The chiral condensate versus /i/(m7r/2) (upper) and versus (lower). 
The dashed curves in the lower figure are drawn to guide the eye and in the lower figure 
they represent the mean field result (62) (Data points taJcen from ref. [84]). 



and of MiTi^ + M2S. The next-to-leading order chiral Lagrangian can be 
written as 

9 

LW=5^L,0i. (63) 

At next-to- leading order, that is p^. one has to take into account the 
one-loop diagrams from the leading-order chiral Lagrangian (52), as well 
as the tree diagrams from the next- to-leading order Lagrangian (63). In 
this perturbative scheme, three Fcynman diagrams contribute to the free 
energy at next- to-leading order (see fig. 7). The one- loop diagram is diver- 
gent in four dimensions. The theory can be renormalized by introducing 
renormalized coupling constants 



A 



d-4 



(64) 



(47r)<*/2 ' 

where A is the renormalization scale and 7i are numbers that can depend 
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Figure 7. Feynman diagrams that enter into the free energy at next-to-leading order. 
The dot denotes the contribution from L (52), and the boxed 4 the contribution from 
(63). The order in the momentum expansion is also given under ea<;h diagram. 



Nj [2,15]. The renormalization can be carried out order by order in the 
perturbation theory. It does not depend on the chemical potential [15]. 

The main technical difficulty at next-to-leading order comes from the 
computation of the one-loop diagram in the Bose condensed phase: Some 
modes are mixed [ 12,15]. Because of this mixing, one-loop integrals may 
be quite complicated. However, we notice that the angle a that appears 
in (58) can be used as an order parameter of the Bose condensed phase. 
Since we want to study the free energy of that phase near the mean-field 
critical eliemical potential jic = -M/2, it is sufficient to compute the one- 
loop integrals for small a, and fj, close to M/2. The free energy is then 
given by 

^^^~cst- (^a2 + (2 + a3)(^-i))a2+ Q-a4)a^ + ... (65) 

The coefficients ai come from the next-to- leading order corrections. They 
are numbers that can be expressed in terms of the renormalized coupling 
constants (64). Their general form is given by 

\k=0 fc=0 / 

They do not depend on the renormalization scale A and can be evaluated 
from the which can in principle be obtained from lattice simulations. 
They are expected to be small (of the order of 0.05 in 3-color QCD [2]). 

The free energy (65) can be analyzed in the same way as a Landau- 
Ginzburg model. The coefficient of is positive. Therefore, there is a 
second order phase transition when the coefficient of vanishes. We thus 
find that the critical chemical potential at next-to-leading order is given by 

Mc = ^M(l-a2) = im^^°, (67) 

where m^^^ is the mass of the Goldstones at next-to-leading order and at 
zero /i. It is remarkable that the next-to-leading order shift in corre- 
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sponds exactly to the next-to- leading order correction of m^. At next- to- 
leading order, we therefore find a second order phase transition at half the 
mass of the lightest particle that carries a nonzero baryon charge. 

We have also calculated the critical exponents at next-to-leading order 
in chiral perturbation theory and find that they are still given by their 
mean-field values [15]. Since = 4 is the critical dimension beyond which 
mean field exponents become valid, this is not entirely surprising. From 
the form of the propagators, we conjecture that the critical exponents are 
given by mean-field theory at any (finite) order in perturbation theory. 

In summary, we find that the next-to-lcading order corrections arc only 
marginal. The main picture obtained from the mean-field analysis is still 
valid at next-to-leading order: A second order phase transition at /Xc = 
m,r/2 with mean-field critical exponents separates the normal phase from 
a Bose condensed phase. 

4.3. Nonzero Temperature 

At the one-loop level in chiral perturbation theory, the influence of the tem- 
perature on the second-order phase transition can also be studied [ 16,88]. 
In order to study the phase transition at nonzero T and /x, we compute 
the free energy of the Bose condensed phase close the critical chemical po- 
tential lie = TO7r/2 at r = 0. The temperature dependence of the free 
energy is solely contained in the 1-loop diagram in fig. 7. Since we are only 
interested in the behavior of the free energy of the Bose condensed phase 
near the phase transition, it is sufficient to compute it for small a, small 
T, and ^ close to 771^^/2. This procedure again leads to a free energy that 
can be analyzed as a usual Landau-Ginzburg model. The minimum of the 
free energy is given by 



where Cj are coefficients that can be computed exactly. For instance, we 



get that C2 = {32V2^F^ - C{l)Vrrh/r^)^ for Nf = 2[ 16]. We find that 

at 




(68) 





(69) 
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both C2 and C4 vanish. Therefore, we find that the second-order phase 
transition fine given by C2 = 0, that is 

Msec(T) = ^ + ^^^C(3/2), (70) 

ends at /xtri- For a larger chemical potential, the phase transition is of first 
order. The second-order phase transition line (70) is in complete agree- 
ment with the semi-classical analysis of a dilute Bose gas in the canonical 
ensemble. This phase diagram has been confirmed by lattice simulations 
[82,85]. 

5. Conclusions 

Below Aqcd the QCD Dirac spectrum both at zero and at a suflaciently 
small nonzero chemical potential is described completely by a suitable chi- 
ral Lagrangian. Below the Thouless energy, i.e. the scale for which the 
Compton wavelength of the Goldstone bosons is equal to the size of the 
box, the Dirac spectrum can be obtained from the zero momentum part of 
this theory. This matrix integral can also be derived from a chiral Random 
Matrix Theory with the global symmetries of the QCD partition function. 
Therefore, below the Thouless energy, the correlations of QCD Dirac eigen- 
values arc given by chiral Random Matrix Theory. At nonzero chemical 
potential the Dirac eigenvalues are located inside a strip in the complex 
plane. Going inside this strip the chiral condensate rotates into a super- 
fiuid Bose-Einstein condensate. A very similar phase transition is found 
for any system with a chemical potential with the quantum numbers of the 
Goldstone bosons. We have discussed in detail the phase diagram for QCD 
with two colors. Because of the Pauli-Giirsey symmetry diquarks appear 
as Goldstone bosons in this theory. We have analyzed the phase diagram 
of this theory to one-loop order and have found a tricritical point in the 
chemical potential-temperature plane. All our results are in agreement with 
recent lattice QCD simulations. 
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